Effect of dopant atoms on local superexchange in cuprate superconductors: A 

perturbative treatment 
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Recent scanning tunneling spectroscopy experiments have provided evidence that dopant impu- 
rities in high-Tc superconductors can strongly modify the electronic structure of the Cu02 planes 
nearby, and possibly influence the pairing. To investigate this connection, we calculate the local 
magnetic superexchange J between Cu ions in the presence of dopants within the framework of the 
three-band Hubbard model, up to fifth order in perturbation theory. We demonstrate that the sign 
of the change in J depends on the relative dopant-induced spatial variation of the atomic levels in 
the Cu02 plane, contrary to results obtained within the one-band Hubbard model. We discuss some 
realistic cases and their relevance for theories of the pairing mechanism in the cuprates. 
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I. INTRODUCTION 

Scanning tunnelling spectroscopy (STS) experiments 
on surfaces of several high-Tc materials [H have discovered 
a host of fascinating phenomena, including checkerboard 
local density of states (LDOS) modulations and inhomo- 
geneous superconductivity with enormous gap modula- 
tions taking place at the nanoscalejj H, 0j 0|- Recently, 
the size of the local gap was found to be positively cor- 
related with simultaneously imaged atomic scale defects, 
thought to be interstitial oxygen dopants Q. This is a 
surprising result, since it had been expected that an oxy- 
gen, which donates two holes to the Cu02 plane, would 
overdope the system and lead to a smaller gap nearby. 
The positive correlation between the positions of the oxy- 
gens and the gap led Nunner et aZ.[3l to suggest that 
the dopants might be increasing the pair interaction lo- 
cally. This could occur if the local electronic structure 
were altered significantly, so as to modify a spin fluctu- 
ation exchange effective interaction, or possibly a local 
electron-phonon coupling constant. For example, it has 
been observed that a strong correlation exists between 
the distance of the apical oxygen from the Cu02 plane in 
high-Tc materials and the critical temperature d, Q, and 
it might be imagined that a modulation of this displace- 
ment by dopant atoms could change the pairing interac- 
tion locally. 

Nunner et al.^ did not assume any specific micro- 
scopic model, but pointed out simply that the general 
assumption of dopants modulating the pair interaction 
could explain a remarkable number of experimental re- 
sults and correlations. Within a generalized inhomo- 
geneous Bardeen-Cooper-ScliriefFer (BCS) pairing model 
adopted in this work, it correctly reproduces the anti- 
correlation of coherence peak height and position, the 
correlation of dopant position with gap size, and the 
detailed frequency dependence of the 0:LDOS(w) cor- 
relation. The theory is still controversial; in almost 



all treatments of disorder in superconductors, impurities 
are assumed to simply scatter electrons as a screened 
Coulomb potential, rather than modulate the pair inter- 
action. There are, however, well-known exceptions [lo|. 
and it is certainly reasonable to expect modulation of 
the pair interactions to be largest in systems like the 
cuprates where the coherence length is small. Strong 
correlations have not been included systematically in the 
theory, although some first steps have been made when J. 
X. Zhu showed explicitly using an inhomogeneous slave- 
boson approach that the proposal of Nunner et al. that 
impurities might modulate the local exchange J was con- 
sistent with the STM observations. 11 1. 



This scenario was investigated beyond the framework 
of mean field theory by Maska et al.\VM. By assuming 
that (i) the cuprate superconductors can be described by 
the t~J model, with the exchange interaction as the main 
pairing mechanism 1^, 14 1 , and (ii) that the presence of 
the dopant atoms induces a position-dependent shift of 
the atomic levels in the Cu02 plane, these authors calcu- 
lated the effective superexchange interaction J between 
copper ions in the presence of dopants from a perturba- 
tion expansion of the one-band Hubbard Hamiltonian up 
to second order. They showed that the diagonal disor- 
der in the plane always leads to an enhancement of J; 
accordingly, with the assumption that pairing is due to 
superexchange, the superconducting gap increases in the 
vicinity of the dopant atoms, in agreement with Ref. 0. If 
true, this conclusion would provide an important, appar- 
ently robust way to connect local atomic displacements 
with the increase of the pairing there, using the results 
of STS. However, recent results by Johnston et al. 
based on cluster model calculations for the three-band 
Hubbard Hamiltonian [l^, which account explicitly for 
the Cu-0 hopping processes, showed instead that elec- 
trostatic modifications due to the presence of the oxy- 
gen dopant locally suppress J. They showed, in addi- 
tion, that electronic coupling to local phonon modes was 
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strongly modified by the dopant, and could enhance J. 
This is consistent with the fact that the gap inhomo- 
geneities are strongly (anti)correlated to a local bosonic 
mode frequency identified in the tunnelling conductance 
by Lee et al.[17|, but still in apparent contradiction to 
the result of Maska et al. 

In view of the present controversy, we investigate here 
the possibility that the Maska et al. result is an artifact 
of an oversimplified model of the electronic structure of 
the Cu02 plane, and analyze the effect of a dopant im- 
purity on J by performing a perturbation expansion on 
the three-band Hubbard model. As in Ref. 13, we as- 
sume initially that the primary effect of the dopant is 
the shift of the atomic energy levels in the Cu02 plane, 
but account for the shifts in O levels as well as Cu. Our 
calculations show that the fifth-order contribution is as 
important as the fourth-order one, in agreement with 
the results for the pure case We find that the sign 
of these contributions is very susceptible to the relative 
dopant-induced spatial variation of the atomic levels; in 
contrast to the single-band case, it may be either posi- 
tive or negative. Finally, we show how the discrepancies 
between Refs. [l^ andll5| may be understood in terms of 
limiting considerations. 



II. MODEL 

Our starting model is the three-band Hubbard Hamil- 
tonian i/nub on the Cu02 plane The three bands 

arise from the hybridization of: the Cu id^2_y2 orbital 
and the two degenerate O 2p orbitals, O 2px and O 2py. 
In the hole representation, -ffnub can be written as: 

i?Hub = ^{ed + V,)d\j,^„+^{ed + l^ + 5i)pl^Pi^a 

i,(T I, a 

<i,l>cr 
<l,n>a 



In Eq. ([T]), d] ^ {di,a) creates (annihilates) a hole with 
spin a in the 3dx2-y2 orbital of a Cu atom at site i. 
Correspondingly, pj ^ {pi.a) creates (annihilates) a hole 
with spin a in one of the two O 2p orbitals at site /. 
Ed is the on-site energy of the Cu 3dj^2_,j,2 orbital, while 
A is the difference between the Cu idj.2_y2 and the O 
2p energies in the pure system, t^^ and describe the 
nearest-neighbor Cu-0 and 0-0 hoppings, respectively. 
Only hoppings within the Cu02 plane are considered. 
The sign of t^^ and t'j^ depends on the relative phase of 
the overlapping 3d^2_y2 and 2p orbitals. Ud (Up) is the 
on-site Coulomb repulsion for a pair of holes on a Cu 
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FIG. 1: (Color online) Energy level diagram for a CU2O5 
cluster illustrating the notation used for the dopant-induced 
shifts of Cu and O atomic energy levels. Each schematic en- 
ergy level is located directly beneath its corresponding atom 
and accordingly colored [black for Cu and cyan (gray) for O 
levels! . 



(O) atom. The presence of a dopant shifts the atomic 
Cu and O energy levels in its neighborhood. Wc denote 
the energy shift for a Cu at position i as Vi and for an 
O at position I between Cu ions at positions i and j as 
Si = 5ij (see Fig. [1]). Besides these shifts, the dopant is 
expected to cause local lattice distortions, which lead to 
the modification of the hopping integrals and t'j^. In 



the present work we neglect this effect, as in Ref. [12 
and concentrate on the effects due to the dopant-induced 
spatial variation of Cu and O atomic energy levels. 



III. PERTURBATION EXPANSION 

The fourth- and fifth-order expressions for the superex- 
change interaction (j'^' and J^^\ respectively) in homo- 
geneous cuprates were derived by Eskes and Jefferson [I^ 
using Rayleigh-Schrodinger perturbation theory: 
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where we have set to zero the Coulomb repulsion Upd 
between a hole on a Cu ion and a hole on the neighboring 
O ion, which we neglect in our calculations for the sake 
of simplicity. 

We consider now the three-band Hubbard model for 
the case with an impurity [Eq. ([T])] in the regime where 
Vi,5i < A,Ud,Up and apply Rayleigh-Schrodinger per- 
turbation theory. We treat the hopping terms in Eq. ([T|) 
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as a perturbation Hi , 

<i,l>a 



<l,n>a' 



(5) 



The ground state of the unperturbed Hamiltonian Hq 
(^Hub = ^^0 + ^^i) corresponds m this case to aU Cu 
atoms occupied by one hole each. This state is 2^-fold- 
degenerate due to the various possible electron spin dis- 
tributions 
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FIG. 2: (Color online) Graphs describing the (a) fourth- and 
the (b) fifth-order hole hopping processes that result in the 
exchange of spins between two Cu atoms. Black [cyan (gray)] 
circles represent Cu [O] atoms. Arrows denote hopping pro- 
cesses, with the accompanying number indicating the order, 
in which the hoppings occur. Symbols a or a stand for the 
spin of the hole. 



where tri, . . . , cttv =T or | and i runs over Cu sites. 

The effective Hamiltonian HeS is calculated as a per- 
turbation expansion in powers For the 
set of states (O, we can ignore many terms of the per- 
turbation series by making use of the fact that the terms 
containing PHiP, where the operator P projects on the 
ground state manifold Eq. will all vanish since it is 
not possible to connect any two states out of the ground 
state manifold [Eq. ^] by a single hopping process. This 
observation leads to the following expression for Hos- 

H^g = EqP + PHiRHiP + PHiRHiRHiP 
+ PHiRHiRHiRHiP 

- ^PHiR^Hi PHi RHi P - i PHi RHi PHiR^HiP 
+ PHiRHiRHiRHiRHiP 

- -PHiRHiR^HiPHiRHiP 
2 

- -PHiRHiPHiR^HiRHiP 
2 

- -PHiR^HiRHiPHiRHiP 
2 

- -PHiRHiPHiRHiR^HiP 
2 

- -PHiR^HiPHiRHiRHiP 
2 

- -PHiRHiRHiPHiR'^HiP, (7) 

where i? = (1 — P)/{Eo — Hq) so that, for a state 
|(A)^{|tTi---aN)}, 
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Eq — Ea, 



(8) 



Eq is the ground state energy of Hq and E^ = ((/i|i7o|'/')- 
Among the terms in i?offi Eq. ([7]), we need only to 
consider those terms that are of the form 



dj „dl =c?, n-di 



<i,j>,<T 



(9) 



with a ~ —a, since the corresponding prefactor deter- 
mines J. The terms of interest result from calculating 
the fourth-order term PHiRHiRHiRHiP and the fifth- 
order term P77ii?ifii?iJii?iJii?iJiP in Eq. dZl). Another 
terms will only add a constant energy term to the effec- 
tive Hamiltonian. 



It is convenient to use graphs for deriving expres- 
sions for J^'^'> and J^^\ One has to consider all pos- 
sible fourth- and fifth-order hopping processes result- 
ing in the exchange of spins between two Cu atoms 
and sum up the corresponding PHiRHiRHiRHiP and 
PHiRHiRHiRHiRHiP expressions. 



There are in total 12 fourth-order graphs (each of the 6 
topologically distinct graphs has two versions that differ 
by flipped spins), one of which is shown in Fig. [2ja). 
In the pure case they reduce to only two terms in J*^^^ 
[Eq. ([3])]. In the case of an impurity (dopant-induced 
spatial variation of Cu and O levels), the terms in the 
sum for J^^^ corresponding to hoppings that start from 
the Cu ion at site i will differ from those corresponding 
to hoppings that start from a Cu ion at site j due to 
the different dopant-induced shifts Vi and Vj. With the 
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is easy to check that the correction terms rj^ and 77, 
vanish when the impurity-induced potentials Vi vanish. 

We note that the sign of the total fourth-order correc- 
tion due to the presence of the dopant, rjij = r{j^^ + ^ij^ 
depends on the sign and the magnitude of Vi and Vj, i.e., 
the actual energy separation between the dopant-shifted 
Cu and O levels and, in particular, for Vi,Vj « A, 
rjij is proportional to Vi + Vj. In general, this result is 
shown diagrammatically in Fig. [3] for a typical set of 
model parameters in the Bi superconductors \Ud = 8.8 




FIG. 3: (Color online) sign(77ij ) diagram for a typical set of the 
model parameters [Ud = 8.8 eV, Up = 4.1 eV and A = 2.92 
eV (Ref. [ish ] in the space of abscissa Vi = Vi — Sij and ordi- 
nate Vj ~ Vj — Sij. White and cyan (gray) regions indicate 
negative and positive total fourth-order correction to J, re- 
spectively. Points (a)-(d) denote the values of Vi and Vj used 
for generating diagrams (a)-(d) in Fig. [4] 
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The sign of this term depends on 



V, 
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eV, Up = 4.1 eV and A = 2.92 eV (Ref. [15^] . In the ^^^^ ^^^^^ ^ei = ^ ,c vf = vf 

space of iii, ■y^, negative and positive 77^ contributions to 



J are shown as white and gray (cyan) areas, respectively. 
For the parameters considered in the cluster calculation 
of Ref. [iSl, the contribution of ?/y is negative (black dot 
(d) in Fig. [S]) and therefore the fourth-order correction 
suppresses J in that case. 

The number of graphs contributing to the fifth-order 
correction to the supercxchange, J^^\ is 120. They all 
involve hoppings to one of the four corner O ions [Fig. [2] 
(b)] so that j(^) will also depend on the corner O atoms 
energy level shifts (5f/, 5^] , (5f , and (5f/ (Fig. [J). For 
example, the correction to the first term of J^^^ in Eq. ^ 



— _ — _ — A*^^ 

(alternatively, = Vj — SfJ, etc. could be considered). 
These six parameters define the sign of the total fifth- 
order correction as well. 

In order to quantify the effect that hopping to the cor- 
ner O atoms has on the superexchange J, we consider 

J - -3 - "3 (siich 
a symmetry is realized when the dopant atom is located 
on the line connecting two Cu atoms [11]). For given Vi 
and Vj, it is then possible to draw a phase diagram of 
the sign of the total correction, /.( ^ + rjij , in the space of 
vl and In Fig. 31 we present, as an example, four 
such diagrams corresponding to different sets of Vi and 
Vj (four points in Fig. [3]). For calculating these diagrams 
we chose tpd=l-2 eV and <pp=0.5 eV as also considered in 
the cluster calculations Ref. Il5|. The local Cu and O site 
energies calculated by Johnston et al. correspond to the 
choice of Wi=-0.13 and fj=0.08 eV in Fig. |l](d). It can be 
concluded from examining the diagrams in Fig. [H that 
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FIG. 4: (Color online) sign(r7,j + diagrams in the space 
of abscissa v'{ and ordinate v^. White and cyan (gray) cor- 
respond to negative and positive values of riij + fiij, respec- 
tively. The point <=-0.23 eV, Wj==-0.05 eV in diagram (d) 
corresponds to the energy levels distribution shown in Fig. 2 
of Ref. m. 



the parameters vf^, wf^, Vj^, and Vj"^ have to be slightly 
larger than Vi and Vj to induce the change of sign of the 
correction to J as compared to the fourth-order result. 



suppression of J. In Fig. [3l the condition Vj — Vi <^ Aioc 
is fulfilled in the vicinity of the Vi = Vj line (on the line, 
Vj — Vi^ 0) and indeed J is increased in the first quarter 
and reduced in the third quarter of the diagram. In the 
second and fourth quarters the relative variation of levels 
of Cu atoms, Vi ~Vj, becomes equally important. The 
one-band model excludes completely the O atoms, thus 
ignoring one of the two microscopic factors (change in the 
Cu-0 energy levels separation and the relative shift of the 
energy levels of two interacting Cu ions) that govern the 
variation of the local superexchange coupling J. 

We would like to note that, as shown by Eskes and Jef- 
ferson [3l for the homogeneous case, even the fifth-order 
perturbation expansion for J is insufficient for quantita- 
tive estimates of J and gives overestimated values com- 
pared with the experimental (and cluster-model calcu- 
lated) values of J. Such trends, naturally, are also to 
be expected in the disordered case. We calculate the 
value of the forth- and fifth-order superexchange cou- 
pling corrections with the same model parameters as used 
in the cluster-model calculations [ll] [Fig. [3] and [4] (d)] 
and find that the sign and the order of magnitude of 
the correction within our calculation, yielding a suppres- 
sion of J by 0(5%), are in good agreement with cluster 
calculations when no modulation in the hopping inte- 
grals is considered, as it is our case here. Consideration 
of other sets of model parameters 22| lead to the same 
relative correction values. 



IV. ANALYSIS AND DISCUSSION 

We are now concerned with trying to answer the follow- 
ing questions. Is the sign of the correction to J caused 
by the impurity uniformly positive, as occurred in the 
one-band calculation of Ref. For physically reason- 
able assumptions regarding the magnitude and spatial 
dependence of the impurity potential for a dopant sit- 
ting several A from the Cu02 plane, can the modulation 
of J be significant at all? 

Regarding the first point, there is a simple argument 
that explains why the superexchange corrections due to 
doping derived from the three-band model can assume 
both positive and negative values while from the one- 
band model one finds that J is always enhanced. Let 
us consider what happens between neighboring Cu and 
O ions when their energy levels shift due to a dopant 
by V and 5, respectively. The local separation between 
their energy levels, which we denoted as A for the ho- 
mogeneous case, varies as Aioc = A — {V — 6). For 
{V — 5) > 0, Aloe decreases compared with A and vice 
versa. Expressing J in terms of Aioc instead of A, one 
sees that for Vj — Vi ^ Aioc [in this limit A in Eq. ([3]) 
and Eq. ([4] can be replaced by Aioc) the variation of lo- 
cal Cu and O levels separation Aioc defines the change 
of J: since Aioc is in the denominator, {V — S) > leads 
to the enhancement of J and {V — 6) < leads to the 



V. CONCLUSIONS 

The question of the impact of a dopant atom on the 
local electronic properties in the Cu02 plane of the 
cuprates has been highlighted by STM measurements, in- 
dicating that dopants correlate with regions of large gap 
0, and the theoretical proposal that the dopant itself 
is enhancing the pairing interaction Q. The appealing 
argument of Maska et al. based on a single-band 

analysis, that the perturbation provided by the dopant 
necessarily enhances the superexchange locally and may 
therefore enhance pairing, has been shown to be a special 
result restricted to one-band systems. Within a three- 
band Hubbard model appropriate to the Cu02 plane, we 
have in this work performed a perturbative calculation to 
fifth order in the hoppings tpd and tpp, and shown that 
the sign of the correction to J can be positive or negative 
depending on the potentials on nearby sites induced by 
the dopant impurity. The typical modulation is of order 
dA/A times the exchange for the homogeneous system, 
where dA is a typical dopant-dependent modulation of 
the local charge transfer energy between Cu and O, and 
A is the homogeneous value of this difference. Using val- 
ues of these shifts obtained from cluster calculations [3| , 
we find that a typical modulation due to an O dopant in 
the Bi-2212 system imaged by STM is of order 5% of the 



6 



homogeneous value. 

It is possible that a more accurate microscopic calcu- 
lation, accounting for the modulations of the hoppings 
and the apical oxygen degrees of freedom neglected here 
may produce a reliable description of this modulation. 
Until then, we have shown that the size and sign of this 
modulation is not universal but depends on details of the 
impurity and local electronic structure. 
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